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ABSTRACT 

A steady, three-dimensional viscous "average pas- 
sage" computer code is used to analyze the flow through 
a compact radial turbine rotor. The code models the 
flow as spatially periodic from blade passage to blade 
passage. Results from the code using varying computa- 
tional models are compared with each other and with 
experimental data. These results include blade surface 
velocities and pressures, exit vorticity and entropy 
contour plots, shroud pressures, and spanwlse exit 
total temperature, total pressure, and swirl distribu- 
tions. The three computational models used are invis- 
cld, viscous with no blade clearance, and viscous with 
blade clearance. It is found that modeling viscous 
effects Improves correlation with experimental data, 
while modeling hub and tip clearances further improves 
some comparisons. Experimental results such as a local 
maximum of exit swirl, reduced exit total pressures at 
the walls, and exit total temperature magnitudes are 
explained by interpretation of the flow physics and 
computed secondary flows. Trends in the computed blade 
loading diagrams are similarly explained. 

INTRODUCTION 

The use of radial turbines is advantageous in 
small propulsion engines due to high efficiencies and 
the large pressure ratios achievable in one stage. 
Furthermore, the use of only one radial turbine stage 
in place of several axial turbine stages can improve 
reliability through reduced part count. However, the 
size and weight of present radial turbine rotor designs 
makes them undesirable for future advanced aircraft 
applications. In addition, low cycle fatigue life is 
a major concern for man-rated applications based on 
present radial rotor sizes and anticipated future tip 
speeds. A compact radial turbine which is 40 to 
50 percent shorter in axial length, 4 percent smaller 
In diameter, and 20 to 30 percent lighter than current 
s tate-of-the-art radial turbines has been designed and 
fabricated by Pratt and Whitney and tested at NASA 
Lewi s Research Center . 


In an effort to understand the complex flow phys- 
ics associated with the compact radial turbine rotor 
operation and to help explain experimental data, a com- 
putational study of the rotor flow field was under- 
taken. A fully three-dimensional code is desired 
because of the highly complex flow associated with 
radial turbine flow. This is especially true for the 
compact radial turbine owing to the short axial dis- 
tance over which the flow is turned. Previous three- 
dimensional computational study of radial turbines is 
sparse, but Choo and Clvinskas (Ref. 1) reported a 
three-dimensional inviscld analysis of a radial turbine 
flow using the DENTON code (Ref. 2). Viscous effects 
such as boundary layer separation and wake development 
were cited as possible causes of predicted underturning 
relative to experimental exit flow angles. Zangeneh- 
Kazemi et al . (Ref. 3) conducted a three-dimensional 
viscous analysis of a radial turbine and compared the 
results with experimental data. It was recommended 
that a tip clearance model be Incorporated for accurate 
simulation of radial turbine flow fields. Good agree- 
ment was found between experimental and computed 
results, especial ly for the calculation accounting for 
tip clearance. However, it is difficult to Interpret 
the comparisons from a turbine design standpoint due to 
the lack of axlsymmetrlc average comparisons. In addi- 
tion, no Inviscid results are presented for comparison 
to the more accurate models. The present study aims to 
produce comparisons with experimental data by consider- 
ing axlsymmetrlc average quantities. Comparisons are 
also made in the present study between three solutions 
of Increasing physical realism: Inviscid, viscous with 

no blade clearance, and viscous with blade clearance. 
Such comparisons are needed to ultimately extend the 
use of three-dimensional viscous turbomachinery codes 
to the turbine design field, because in this way the 
inaccuracies of the Inviscid and nonclearance assump- 
tions can be Isolated and analyzed. 

In the present study, a three-dimensional viscous 
computer code of Adamczyk (Refs. 4 to 6) is used. This 
code Is an "average passage” code which models the flow 
as spatially periodic from blade passage to blade pas- 
sage in a particular blade row and can be used to model 
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multiple blade rows. This Is particularly advantageous 
in multistage turbomachi nery calculations, where this 
averaging process keeps computational time from becom- 
ing prohibitive. For a single blade row, the code is 
essentially identical to a fully three-dimensional vis- 
cous code. A tip clearance model is included in the 
analysis, as recommended by Zangeneh-Kazemi et al . In 
addition, a blade hub clearance model is included which 
models the backface clearance typically found in radial 
turbine rotors. 

ROTOR DESCRIPTION 

The compact radial turbine rotor Is shown in 
Fig. 1, The rotor has 14 solid blades attached to a 
solid disk. The rotor has radial inflow and axial out- 
flow and is scalloped; that is, the disk does not 
extend to the leading edge of the blades. Instead, the 
flow is bounded by the stationary backface in this 
region. Likewise, the stationary aft centerbody is in 
contact with the flow downstream of the rotor. Figure 2 
details the relevant rotor features. The shroud clear- 
ance Is about 0.046 cm (0.018 In.) at the leading edge 
and 0,020 cm (0.008 in.) at the trailing edge. The 
backface clearance is about 0.038 cm (0.015 in.). The 
rotor diameter is 0.367 m (14.5 in.) to the blade lead- 
ing edge. 

Figure 3 illustrates the 40 percent reduction In 
axial length achieved with the compact radial turbine 
design compared to a state-of-the-art radial turbine of 
the same class. The reduction in length holds several 
advantages. In addition to the size and weight reduc- 
tions important for aircraft applications, a shorter 
rotor reduces the unsupported shaft length between 
bearings, resulting in improved rotor dynamics. 

The experimental study of the compact radial tur- 
bine was conducted in the NASA Lewis Warm Turbine Test 
Facility under a joint program with Pratt and Whitney. 
Due to the confidential nature of the program, this 
data has not been fully published and disclosure of 
turbine efficiencies is restricted. Preliminary test- 
ing consisted of a tare test on the basic rig with a 
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Baseline Rotor 


Figure 1 - Compact radial rotor 


Figure 3 - Compact rotor profile comparison 

bladeless rotor. This was done to measure the losses 
associated with the rig Itself; namely windage, bear- 
ing, and seal losses. Next, the rotor was Installed 
into the rig, and testing was conducted over a matrix 
of test conditions. This matrix consisted of a range 
of pressure ratios and rotor speeds centered about the 
design point. The design point for the turbine was 
20 000 rpm, 5.0 pressure ratio, 2.86 kg/s (6.30 Ib/s) 
flow rate, and 478 K (860 °R) rotor total Inlet temper- 
ature. This corresponds to a Reynolds number of 
1.90xl0 7 based on rotor diameter and inlet conditions. 
Rotor exit surveys were conducted, which measured flow 
angles, total temperatures, and total pressures versus 
percent span. In addition, shroud static pressures 
were recorded versus meridional distance. All data was 
recorded at a steady-state condition. The design point 
data from these experiments is used for comparison to 
calculation in this report. 
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CODE DESCRIPTION 

The code used to analyze the compact radial tur- 
bine is a three-dimensional viscous "average passage" 
code (Refs. 4 to 6). The code uses three averaging 
processes to obtain the average passage equations from 
the full Navier-Stokes equations. These are ensemble 
averaging, time averaging, and passage-to-passage aver- 
aging. Ensemble averaging yields the familiar 
Reynolds-averaged Navier-Stokes equations. This aver- 
aging process removes the random fluctuations from the 
flow field. Time averaging further simplifies the flow 
equations by removing the deterministic unsteady flow 
component which occurs because of the relative movement 
of different blade rows. Finally, the equations are 
passage-to-passage averaged, which removes the effect 
of different blade count of stator and/or rotor blade 
rows. In this procedure, the flow is assumed to be 
spatially periodic from blade passage to blade passage 
in a particular blade row. For multistage turbomachin- 
ery, with stator blade rows and/or rotor blade rows 
having different numbers of blades, this will not gen- 
erally be true. The average passage equations obtained 
by the three averaging processes are discretized using 
a control volume approach in a cylindrical coordinate 
system. Using this control volume approach, all flow 
quantities are cel 1 -centered . To obtain values of flow 
variables at cell surfaces, an averaging process is 
used which results in second-order accuracy. In the 
viscous solutions, a Bal dwi n-Lomax turbulence model Is 
employed. The turbulence model is updated after every 
10 iterations. 

The discretized average passage equations obtained 
by the control volume approach are solved using a four- 
stage Runge-Kutta scheme. Second and fourth difference 
smoothing operators are used to avoid alternate point 
decoupling in the solution. The Riemann Invariants C + 
and C~ are used to calculate the radial velocity compo- 
nent at the upstream boundary. Constant entropy Is 
assumed from an upstream state to the inlet boundary, 
and this is used to determine the inflow temperature, 
pressure, and density. At the downstream boundary, the 
static pressure at the hub is specified and radial 
equilibrium Is enforced across the span. All other 
downstream flow quantities are extrapolated from the 
interior of the flow domain. At the blade, hub, and 
shroud surfaces, an extrapolated normal pressure gradi- 
ent boundary condition is employed in conjunction with 
an adiabatic wall. The hub velocity boundary condition 
was adapted for the scalloped rotor configuration. In 
the viscous solutions, the computational cells which 
lie in the stationary backface and aft centerbody 
regions employ a zero velocity boundary condition. All 
other hub cells use the rotor speed as the velocity 
boundary condition, which is typical. The shroud is of 
course stationary in the absolute frame of reference 
and is so modeled. For the inviscld solution, the 
absence of viscous forces renders rotation of axi sym- 
metric surfaces Inconsequential. The mesh boundaries 
extending from the blade leading and trailing edges to 
the inflow and outflow boundaries, respectively, are 
periodic boundaries. At these boundaries, periodicity 
1$ enforced by setting the flow variables In the compu- 
tational cell outside the computational domain equal to 
those at the corresponding periodic location. 

The code contains a one-dimensional inviscld flow 
model for the clearance gap flow. This model conserves 
mass, momentum, and total enthalpy across the clearance 
gap. As a result, the pressure difference across the 
gap is zero. This is a rather simplified representa- 
tion of the flow In the clearance region, but it does 
result in a flow across the clearance gap which serves 
to reduce the loading in the tip region. This model is 


simply implemented by enforcing the previously stated 
periodic boundary conditions. This model for the flow 
in the clearance gap was applied to the entire blade 
tip and to the hub backface clearance. 

METHOD OF ANALYSIS 

The average passage code was first modified for use 
with radial turbomachinery. Among these modifications 
was the changing of the upstream velocity boundary con- 
dition from axial to radial flow and the replacement of 
various axial velocity approximations by the true 
velocity. The code was then in an appropriate form to 
analyze the compact radial turbine rotor. Although the 
code Is capable of analyzing the entire turbine stage, 
this was not done as part of this analysis. Instead, 
it was assumed that the flow upstream of the rotor was 
axi symmetric and spanwise constant. In a true three- 
dimensional stage calculation, neither of these assump- 
tions would be made. The experimental inlet flow was 
not spanwise constant, since boundary layers formed on 
the hub and shroud surfaces. However, spanwise con- 
stancy was assumed in the calculation for simplicity. 

A three-dimensional computational grid was gener- 
ated for the rotor blade passage using an interactive 
H-grid generator. The hub and suction surface grid Is 
shown In Fig. 4. The grid size is 109 by 33 by 33 grid 
points, with 109 grid points In the flow direction. 

The grid includes densely-packed grid points near the 
solid boundaries to adequately resolve viscous effects. 
The wall spacing is from 0,0025 to 0.0050 cm (0.0010 to 
0.0020 in.) over most of the blade, hub, and shroud 
surfaces. The grid extends from about 13 percent of 
midspan chord upstream of the blade leading edge to 
about 27 percent of midspan chord downstream of the 
blade trailing edge. The midspan chord Is defined as 
the arc length of the midspan grid line on the blade 
surface. Grid points are packed densely in the stream- 
wise direction around the elliptic leading and trailing 
edges of the blades. This Is necessary to adequately 
resolve the accelerating flow in these regions. The 
grid spacing for the first two grid points in the 
streamwlse direction around the blade leading and 
trailing edges is about 0.011 cm (0.0043 in.). There 
are two grid points In the tip clearance region and 
three in the hub clearance region. Only one grid was 
employed in the calculations, so the effect of differ- 
ent grid densities could not be studied. 

The code requires specification of upstream total 
temperature and pressure, upstream swirl velocity, 
downstream static pressure at the hub, and rotor speed. 
Downstream static pressure, upstream total temperature, 
and rotor speed were taken to match experimental condi- 
tions at the design point. The upstream swirl velocity 
was taken to match design conditions based on the sta- 
tor exit blade angle. The upstream total pressure was 
calculated based on a 2 percent reduction from the 
stator inlet value. The mass flow rate is calculated 
by the code. A fuel-to-alr ratio of 0.0042 was used 
to determine the gas constant of 288.1 mws 2 K 
(1723 ftw$2 °R) . Since constant specific heats were 
assumed in the calculation, a mean value of 1.39 was 
used for the specific heat ratio. 

A two-dimensional inviscld hub-to-shroud midchan- 
nel stream surface solution was obtained for the rotor 
flow field using the MERIDL code (Ref. 7). This solu- 
tion was Interpolated to the three-dimensional grid for 
use as a starting solution for the three-dimensional 
average passage calculation. Three solutions of the 
rotor flow field were made using varying computational 
models. The first solution employed the Inviscld form 
of the code. The second solution employed the viscous 
turbulent form of the code. This solution attempted to 
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Figure 4 • Computational grid 


accurately model the rotor backface and aft centerbody 
by specifying these regions as stationary. The third 
solution differed from the second solution only in that 
it included the effects of the hub and tip clearances* 
The hub clearance Is due to the separation distance 
between the rotor backface and the blades in the radial 
portion of the rotor. All three solutions employed the 
same computational grid. The code was run for each of 
the three cases until the root-mean-square average of 
the residuals at each grid point had decreased by at 
least three orders of magnitude from the starting solu- 
tion. The residual was defined as the time derivative 
of the gas density at a specific location. 

The computations were performed on a CRAY-2 super- 
computer. The inviscid and viscous solutions required 
approxi mate 1 y 2000 and 4000 iterations, respectively, 
to reach a converged solution from a rough starting 
solution. This corresponds to approximately 1.75 and 
4.25 CPU hours of run time, respectively. The computa- 
tional grid, executable code, and flow solution required 
about 2.5 klloblocks of combined memory. 

RESULTS AND DISCUSSION 

Figures 5 to 7 correspond to the viscous nonclear- 
ance solution. Figure 5 presents the relative velocity 
vectors at the grid points adjacent to the hub and 
blade suction surfaces. Over most of the blade suction 
surface, the spanwise component of the flow adjacent to 
the blade is directed toward the shroud. However, near 
the blade tip the flow vectors point away from the 
shroud. This is due to viscous effects at the shroud. 

In the rotating frame of reference, the fluid adjacent 
to the shroud is dragged by the shroud into the suction 
surface of the blade. This causes the flow near the 
shroud to curl into a counterclockwise vortex (as 
viewed from downstream) known as the scraping vortex. 
This scraping vortex Is counter to the direction of the 
passage vortex. The clockwise passage vortex is pro- 
duced by the components of streamline curvature and 
Coriolis acceleration in the direction of absolute vor- 
ticity of the flow and is an inviscid effect. Figure 6 


shows the counterrotating vortices In more detail. 

This Is a contour plot of axial vorticity on the grid 
plane containing the blade trailing edges as viewed 
from a downstream viewpoint. The vorticity Is non- 
dimensionali zed by the ratio of reference velocity to 
reference length. The reference velocity is the stag- 
nation sonic velocity at the rotor Inlet divided by the 
square root of the specific heat ratio and the refer- 
ence length is the rotor diameter. The blade suction 
surface Is to the left of the figure. The vortices are 
adjacent to the suction surface near the shroud. It 
can be seen from the sense of rotation of the vortices 
that the fluid between them is driven away from the 
wall. This Indicates that the boundary layer fluid on 
the suction surface is pulled away from the blade sur- 
face Into the main flow by the vortices. Since the 
boundary layer flow contains high loss fluid, this vor- 
tex action brings higher loss fluid into the main flow. 
This is confirmed by Fig, 7, which is a plot of entropy 
on the same exit plane. The entropy Is nondlmenslonal- 
ized by the gas constant and the zero entropy datum is 
at the rotor inlet. The high entropy (high loss) flow 
is centered between the two counterrotating vortices. 
Figure 6 also shows regions of high vorticity along the 
blade pressure surface. This is associated with the 
origin of the viscous wake. Figure 7 establishes this 
region as high entropy, high loss flow. 

The final converged solution for the viscous blade 
clearance model represents the most physically realis- 
tic model used in the analysis of the rotor. Figures 8 
to 11 correspond to this solution. Figure 8 shows the 
relative velocity vectors at the grid points adjacent 
to the hub and blade suction surfaces. Over most of 
the blade suction surface, the spanwise component of 
the flow adjacent to the blade Is directed toward the 
shroud. This effect Is again due to the passage vortex 
which rotates In the clockwise direction. The vectors 
near the hub and tip which are in the circumferential 
direction represent the leakage flow in the blade 
clearance regions. This leakage flow reduces the effi- 
ciency of the machine since the fluid does no work on 
the blade. 



Flour. 7 - Non-cl.aranc* vl.cou. aolutlon axil entropy Flgur . a . clearance vlacou. aolutlon relative velocity vector. 


Figure 9 shows the axial component of vorticity on destroyed by the clearance flow. However, another 
the grid plane containing the blade trailing edges. mechanism causes loss in the clearance solution which 

The vorticity is nondimensional ized by the ratio of is not present in the nonclearance solution. The high 

reference velocity to reference length. The reference loss tip clearance flow passes over the blade tip along 

velocity Is the stagnation sonic velocity at the rotor the entire length of the blade and penetrates into the 

Inlet divided by the square root of the specific heat flow along the shroud toward the adjacent blade. The 

ratio and the reference length Is the rotor diameter. high loss In the clearance flow Is produced by viscous 

The tip clearance flow is driven by pressure differ- effects at the shroud. In the rotating frame of refer- 
ences from the pressure to the suction surface (left to ence, the fluid very close to the shroud rotates with 

right In Fig. 9) and causes the vorticity to be nega- the shroud due to viscous effects. This causes the 

tlve (clockwise) over most of the region near the sue- relative circumferential velocity very near the shroud 

tlon tip corner. The counterrotating vortex pattern to be In the opposite direction to the clearance flow 

present in the nonclearance viscous solution is The resultant high shear region generates loss in the' 
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Figure 9 - Clearance viscoua solution exit axial vorticlty 


clearance flow as it travels along the shroud. This is 
illustrated in Fig. 10, which is a plot of entropy on 
the same exit plane as Fig. 9. The entropy is nondi- 
mensional ized by the gas constant and the zero entropy 
datum is at the rotor inlet. The high entropy region 
adjacent to the shroud and slightly away from the suc- 
tion surface is due to the high loss clearance flow 
which accumulates In this region. Figure 9 indicates 
regions of large vorticlty near the blade pressure sur- 
face, as was seen in Fig. 6 for the nonclearance vis- 
cous solution. As previously noted, this corresponds 
to the beginning of the turbulent wake, and is 
reflected as Increased entropy In Fig. 10. 

It appears from Figs. 6, 7, 9, and 10 that the 
boundary layer on the blade pressure surface Is thicker 
than that on the blade suction surface. These figures 
are plotted on the grid plane containing the blade 
trailing edges. From Fig. 4 it can be seen that this 
plane is curved, particularly near the blades. This 
curvature causes more of the grid plane to pass through 
the blade aft wake on the pressure side. This is con- 
firmed by Fig. 11, which is a contour plot of relative 
critical velocity ratio on the bl ade-to-bl ade midspan 
surface viewed from an oblique angle. It is the grid 
plane curvature effect which is responsible for the 
illusion of a thicker boundary layer on the blade pres- 
sure surface in Figs. 6, 7, 9, and 10. Figure 11 indi- 
cates that this Is not the case, as the suction surface 
boundary layer is clearly thicker than that of the 
pressure surface. 

It can be seen from Fig. 11 that the flow acceler- 
ates from a relative velocity ratio of about 0.3 at 
inTet to approximately sonic relative velocity at the 
exit. The turbulent wake behind the blade trailing 
edge Is characterized by the large velocity gradient In 
this region. The blade wake direction follows the 
blade exit angle closely. Another region of large 
velocity gradients is on the blade suction surface near 
the leading edge. Here the flow is rapidly expanding 
and accelerating around the leading edge after stagnat- 
ing on the nearby pressure surface. 

A very small separation region was detected in the 
nonclearance viscous solution. This region is located 
at the intersection of the hub and the blade suction 
surface, where the hub first begins to turn from the 
radial direction. However, this separation region was 


Figure 10 - Clearance viscoua solution exit entropy 



Figure 11 - Clearance viscous solution blade-to- 
blade midspan critical velocity ratio 


not encountered in the clearance viscous solution. 

This is due to the blade hub clearance flow which 
entrains the boundary layer flow in the radial portion 
of the hub, and keeps the flow from separating. 

Using the 5.0 pressure ratio condition, the invis- 
cid code predicted a mass flow rate of 2.87 kg/s 
(6.33 Ib/s). The viscous nonclearance code predicted 
a mass flow rate of 2.79 kg/s (6.16 lb/s), and the vis- 
cous clearance code predicted a mass flow rate of 
2.78 kg/s (6.12 Ib/s). The experimental value was 
2.86 kg/s (6.30 lb/s). This indicates that the viscous 
solutions underpredicted the flow rate due to the 
blockage effect of the wall boundary layers. Mass flow 
rate was conserved for each of the three calculations 
to within 0.5 percent from Inlet to exit. Doubling 
the artificial viscosity terms from the typical values 
reduced the converged mass flow rate by less than 
0.5 percent. 
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Figures 12 to 15 present computed quantities In 
comparison to experimental results, which are by their 
nature averages. In order to make this comparison, 
the computed local quantities must be averaged In the 
bl ade-to-bl ade direction. The local density, pressure, 
Internal energy, and three components of momentum were 
circumferentially averaged on a flow area basis along 
blade-to-blade grid lines. These axl symmetric average 
quantities were then used to calculate any derived axl- 
symmetric average quantities such as the total pres- 
sure, total temperature, and flow angle. These are the 
quantities plotted In Figs. 12 to 15. 

Figure 12 compares the rotor exit swirl flow 
angles for the three previously discussed solutions 
with experimental data. The exit flow angle Is defined 
as the angle between the absolute exit velocity and the 
axial direction, neglecting the radial velocity compo- 
nent. The direction of rotor rotation Is considered 
positive. The Invlscld solution fails to predict the 
local maximum at 75 percent span Indicated by the data. 
However, the nonclearance viscous calculation does pre- 
dict this maximum. This can be explained by the coun- 
terrotating vortex behavior present in the viscous 
solution. The fluid between the passage and scraping 
vortices tends to underturn; that Is, Its velocity In 
the rotational direction Is greater than the main flow. 
This corresponds to the maxima at 75 percent span In 
Fig. 12. The clearance viscous solution however, does 
not predict this behavior as well. This Is due to the 
tip clearance flow which dominates the scraping vortex 
and causes the entire region to have negative (clock- 
wise) axial vortldty. Thus, underturning does not 
occur at 75 percent span, but much nearer the shroud. 
This effect causes the two viscous solutions to bracket 
the data at the shroud. Such bracketing of the data 
may indicate that the tip clearance In the axial por- 
tion of the rotor Is overestimated In the clearance 
viscous solution. The tip clearance model considers 
the blade tip to be along a grid line. However, the 
grid Is much denser from hub to tip in the radial por- 
tion of the rotor due to the much smaller span there. 
Thus, it was necessary to underpredict the clearance In 
the radial portion and overpredict the clearance in the 
axial portion to minimize the error caused by the vary- 
ing clearance model. This problem was not encountered 
with the hub clearance, since the hub clearance is 
entirely in the radial region of the rotor. It Is 
unclear why all solutions predict underturning relative 
to the experimental data In the 30 to 40 percent span 
region. Choo and Civinskas (Ref. 1) noted a similar 
trend in their invlscld analysis of a radial turbine 
for span locations less than 40 percent. This was 
explained by viscous effects such as boundary layer 
and wake development which would not be modeled by an 
Inviscid analysis. However, the viscous and Invlscld 
analyses in this study predict a nearly Identical rotor 
exit flow angle in the 30 to 40 percent span range. 

Figure 13 compares the rotor exit total pressures 
for the three solutions with experimental data. The 
pressures are normalized with respect to the rotor 
Inlet total pressure. The viscous solutions correctly 
predict a reduction in total pressure very near the hub 
and shroud. This Is due to the boundary layer flow 
near the walls. The clearance solution predicts a drop 
In total pressure at about 80 to 90 percent span. This 
may again be due to an overestimation of the tip clear- 
ance in the axial section of the rotor. However, some 
loss of total pressure In the shroud region is expected 
due to tip clearance losses. 

Figure 14 compares the rotor exit total tempera- 
tures for the three solutions with experimental data. 

The temperatures are normalized with respect to the 
rotor Inlet total temperature. Due to the vortex 


behavior and boundary layer losses present In the non- 
clearance viscous solution, It predicts a higher exit 
total temperature (higher loss) than the Invlscld solu- 
tion at about 80 percent span. This brings the solu- 
tion closer to the data, but the average exit total 
temperature is still lower than the experiment. This 
indicates an underprediction of loss. When the clear- 
ance model is Introduced, the exit total temperature 
again Increases and corresponds more closely to the 
data. This Implies that the clearance model produces a 
loss relative to the nonclearance viscous calculation. 

Figure 15 compares the axl symmetric average shroud 
static pressures for the three solutions with experi- 
mental data. The pressures are normalized with respect 
to the rotor inlet total pressure. All three solutions 
agree well with each other. This indicates that the 
shroud pressure distribution Is governed more by the 
Inlet and exit boundary conditions and the rotor geome- 
try than by viscous and clearance effects. However, 
the solutions differ from the data in the 0 to 50 per- 
cent meridional distance range. This is perhaps due to 
an inaccurate estimation of rotor pressure ratio. The 
exit pressure is accurately matched, as expected. 

Figures 16 to 18 show the static pressures on the 
blade suction and pressure surfaces at the hub, mid- 
span, and tip, respectively, for the three solutions. 
The hub, midspan, and tip stations are the same for 
each solution. The hub and tip stations were taken as 
the first grid points on the blade surface for the 
clearance solution. This corresponds to the fourth 
grid point from the hub and the third grid point from 
the tip, respectl vely. The pressures are normalized 
with respect to the rotor Inlet total pressure. In the 
radial portion of the rotor (approximately 0 to 30 per- 
cent meridional distance) the blade loading Is fairly 
constant from hub to tip. This Is because the wheel 
speed does not vary from hub to tip In this portion of 
the rotor. However, In the axial portion of the rotor 
(approximately 60 to 100 percent meridional distance) 
the loading increases from hub to tip. This is typical 
axial turbine behavior, and arises due to the Increase 
In wheel speed from hub to tip. 

At the tip, the clearance viscous calculation pre- 
dicts a higher loading than the nonclearance viscous 
calculation. To a lesser extent, this Is also true In 
the hub clearance region. For the nonclearance viscous 
solution, the suction surface of the blade passes over 
the shroud with no leakage. Therefore the blade tip 
becomes a stagnation line with a static pressure equal 
to the relative total pressure. However, when the 
clearance is introduced, the blade tip is no longer a 
stagnation line, since there is leakage. Thus the 
blade suction tip static pressure is lower than the 
total pressure, and lower than the nonclearance value. 

In the Invlscld solution, the shroud rotation has no 
effect, so as in the clearance viscous solution the 
suction surface tip static pressures do not reach the 
total pressure. It is interesting that the inviscid 
solution agrees with the clearance viscous solution 
better than the viscous nonclearance solution. This 
Is due to the physically unrealistic condition which 
exists in the viscous nonclearance solution. It Is 
physically unrealistic to assume that the fluid Is 
bound by the no-slip condition to both the shroud and 
to the blade suction surface, yet flow cannot escape 
between the two surfaces as clearance flow. However, 
this Is a highly localized effect, limited to a region 
very near the clearance. Over most of the blade sur- 
face, the clearance solution predicts a reduced loading 
such as that seen in the midspan plot. Thus, the over- 
all blade loading is smaller for the clearance viscous 
calculation than for the nonclearance viscous calcula- 
tion, as expected. 
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Figure 14 - Rotor exit total temperature versus percent span 
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Figure 15 - Axisymmetrlc shroud pressure versus meridional distance 



NORMALIZED STATIC PRESSURE NORMALIZED STATIC PRE! 



Figure 16 • Blade loading diagram at hub 



Figure 17 - Blade loading diagram at mldapan 
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Figure 18 - Blade loading diagram at tip 


CONCLUDING REMARKS 

The three-dimensional average passage code was 
modified for use with radial turbomachinery. It was 
applied to analyze the flow field in an advanced com- 
pact radial turbine rotor. Progressively Improved com- 
putational models were employed, Including an Inviscid 
solution, a viscous solution without clearance, and a 
viscous solution with clearance. Comparisons of the 
calculated flow fields with test data showed good 
agreement with total pressure and total temperature, 
but rather poorer agreement with flow angle. The 
agreement with rotor exit total temperature was best 
for the most realistic model, namely the viscous blade 
clearance model. Both viscous solutions better pre- 
dicted experimental trends In rotor exit flow angle 
than did the Inviscid solution and better predicted hub 
and shroud effects In the rotor exit total pressure. 
However, the clearance viscous solution deviated from 
the flow angle and total pressure data near the shroud. 
This was perhaps due to an overestimation of tip clear- 
ance in the axial portion of the rotor. Near the mid- 
span region the viscous clearance solution showed the 
poorest agreement with the measured flow angle and 
total pressure. 

The average passage code allows for the Interpre- 
tation of secondary flows. It also brings a better 
Insight Into the flow physics of turbomachinery, which 
is particularly complex In radial machines such as the 
compact radial turbine rotor. This Insight can lead to 
an explanation of experimental trends, which makes the 
code a powerful new analysis tool. With experience, 
the knowledge of secondary flow patterns can lead to 
improved turbomachinery designs and Increased engine 
efficiencies. 
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